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ABSTRACT 

Aims. An asteroseismological study of PG 1159—035, the prototype of the GW Vir variable stars, has been performed 
on the basis of detailed and full PG1159 evolutionary models presented by Miller Bertolami & Althaus (2006). 
Methods. We carried out extensive computations of adiabatic gr-mode pulsation periods on PG1159 evolutionary models 
with stellar masses spanning the range 0.530 to O.741M0. These models were derived from the complete evolution of 
progenitor stars, including the thermally pulsing AGB phase and the born-again episode. We constrained the stellar 
mass of PG 1159—035 by comparing the observed period spacing with the asymptotic period spacing and with the 
average of the computed period spacings. We also employed the individual observed periods reported by Costa et al. 
(2007) to find a representative seismological model for PG 1159—035. 

Results. We derive a stellar mass in the range 0.56 — O.SOM© from the period-spacing data alone. We also find, on 
the basis of a period-fit procedure, an asteroseismological model re prese ntative of PG 1159—035 that reproduces the 
observed period pattern with an average of the period differences of 511^ — 0.64 — 1.03 s, consistent with the expected 
model uncertainties. The model has an effective temperature Tcff — 128 OOOI2 eoo a stellar mass 

A/* = 0.565tHo9Afe> 

a surface gravity logy = 7.42lg j2, a stellar luminosity and radius of log(Lt/L0) = 2.15 ± 0.08 and log(_Rt/_R0) — 
— 1. 62^0 03, and a He-rich envelope thickness of Menv = O.O17M0. The results of the period-fit analysis carried out 
in this work suggest that the surface gravity of PG 1159—035 would be 1 a larger than the spectroscopically inferred 
gravity. For our best-fit model of PG 1159—035, all of the pulsation modes are characterized by positive rates of period 
changes, at odds with the measurements by Costa & Kepler (2007). 

Key words, stars: evolution — stars: interiors — stars: oscillations — stars: variables: other (GW Virginis) — white 
dwarfs 



• ^ ' 1. Introduction 2001, Blocker 2001, Lawlor & MacDonald 2003, Althaus et 

al. 2005, Miller Bertolami et al. 2006 for recent references). 

^ ^ PG 1159-035 (GW Virginis) is the prototype of both the pQ ii59_035 is characterized by T,s = 140 000 ± 5 000 

. .C^ , class of pulsatmg PG1159 stars —commonly known as GW ^ and log^ = 7.0 ± 0.5 (Werner et al. 1991). Since the 

Vir or DOV variables— and the spectral class of hot discovery of its photometric variations by McGraw et al. 

hydrogen-deficient, post-asymptotic giant branch (AGB) ^^979)^ ^-^^ ^^^^ been the focus of intensive scrutiny, 

stars with surface layers rich m helium, carbon and oxygen impressive pulsational study, Winget et al. (1991) 

—the PG1159 stars (see Werner & Herwig 2006). GW Vir (WEA91) identified 122 peaks in the periodogram (Fourier 

variables exhibit multiperiodic luminosity fluctuations with spectrum) of PG 1159-035 obtained with the Whole Earth 

periods m the range 5 - 50 minutes, attributable to non- Telescope (WET; Nather et al. 1990). The peaks were 

radial gravity mode pulsations driven by the K-mechamsm attributed to high order nonradial g-modes with periods 

due to partial ionization of carbon and oxygen m the outer between « 300 - 1000 s, and with period spacings of 

layers (Starrfield et al. 1983; Gautschy et al. 2005; Corsico ^jjO^^ ^ 21.50 ± 0.03 s and AUf^^ = 12.67 ± 0.03 s. 

et al. 2006; Quirion et al. 2007). PG1159 stars are con- ^he mode with the largest amplitude (7.2 mma), which 

sidered the evolutionary hnk between post-AGB stars and identified as a £ = 1, m = +1, has a period of « 516 

most of the hydrogen-deficient white dwarfs. It is accepted ^gj^g ^ larger data set from different years (1983, 

that these stars have their origin m a born-agam episode ^935^ ^ggg^ ^993^ ^^d 2002) and improved data reduction 

induced by a post-AGB hehum thermal pulse (see Herwig ^^d data analysis, Costa et al. (2007) (CEA07) identified 

76 additional pulsation modes, enlarging to 198 the total 



Send offprint requests to: A. H. Corsico number of pulsation periods for PG 1159-035, and placing 

* Member of the Carrera del Investigador Cientffico y if as the star with the largest number of modes detected 
Tecnologico, CONICET, Argentina. after the Sun. CEA07 found AII^^^ = 21.43 ± 0.03 s and 



Fellow of CONICET, Argentina. ^^1=2 = 12.38 ± 0.01 
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Fig. 1. Our PG1159 full evolutionary tracks in the 
logTcft — log 5 plane, with stellar mass (from right to left): 
0.530,0.542,0.565,0.589,0.609,0.664, and 0.742AfQ. The 
location of PG 1159—035 according to spectroscopy is indi- 
cated with a star symbol. The hollow squares (diamonds) 
display the predictions of the asymptotic (computed) pe- 
riod spacings fi j3.1l andi |3.2p . The filled circle shows the lo- 
cation of PG 1159—035 according to the period fitting pro- 
cedure O. Thin dotted lines define three different regimes 
for the theoretical rate of period change: models exhibit 
positive (negative) rates of period change below (above) 
the lower (upper) line, and positive and negative rates (be- 
tween both lines). The hot boundaries of the theoretical 
dipole {£ — 1, solid line) and quadrupole {£ — 2, dashed 
line) GW Vir instability domain — according to Corsico et 
al. (2006) — are also shown. For details, see text. 



The determination of the stellar mass of PG 1159—035 
has been the subject of numerous investigations. The stel- 
lar mass of pulsating pre-white dwarfs can be constrained 
from asteroseismology — the asteroseismological mass — ei- 
ther through the observed period spacing (see, for instance, 
Kawaler & Bradley 1994 (KB94); Corsico & Althaus 2006) 
or by means of the individual observed periods (see, e.g., 
KB94, Corsico & Althaus 2006, Corsico et al. 2007ac). The 
study of WEA91 found a mass of 0.586 ± 0.003Mq for 
PG 1159—035 on the basis of the observed period spac- 
ing for £ = 1 and £ = 2 modes and an asymptotic anal- 
ysis based on the PG1159 structure models of Kawaler 
(1986, 1987, 1988). From a detailed period fit, KB94 found 
M* = O.59±O.O1M0 for PG 1159-035. Corsico & Ahhaus 
(2006) found a stellar mass of 0.56A/o from a period fit- 
ting procedure based on PG 1159 evolutionary models with 
several masses created artificially from the full sequence of 
0.589Mq of Althaus et al. (2005). Finally, with the enlarged 
set of periods for PG 1159-035, CEA07 obtain a value of 
M^, = 0.59 ± O.O2M0 by using the parameterization of the 
asymptotic period spacing of KB94. 

The total mass of PG1159 stars can also be estimated 
through the comparison of the spectroscopic values of Tcs 



and g with evolutionary tracks — the spectroscopic mass. 
On the basis of the evolutionary tracks of O'Brien & 
Kawaler (2000), Dreizler & Heber (1998) derived a stel- 
lar mass of 0.54 ± O.IMq for PG 1159-035. The most re- 
cent determination is that of Miller Bertolami & Althaus 
(2006), who also derived a stellar mass of 0.54 ± O.IMq 
on the basis of PG1159 evolutionary models with different 
stellar masses that take fully into account the evolutionary 
history of their progenitor stars, particularly the thermally 
pulsing and born again phases. 

The discrepancy between the asteroseismological and 
the spectroscopic mass of PG 1159—035 has been par- 
tially alleviated by the employment of Miller Bertolami & 
Althaus (2006) PG1159 evolutionary models and the av- 
erage of the computed period spacings in the determina- 
tion of the asteroseismological mass. In fact, at the effective 
temperature of PG 1159—035, the average of the computed 
£ — 1 period spacings is consistent with a stellar mass of 
« 0.558Mq (Corsico et al. 2006). 

The aim of this paper is to present a detailed asteroseis- 
mological study of PG 1159—035 based on the evolutionary 
models of Miller Bertolami & Althaus (2006) and the ob- 
servational data of CEA07. The present study is the third 
time that such consistent PG1159 evolutionary models are 
employed in individual asteroseismological modeling of pul- 
sating PG1159 stars, the first application being the study 
carried out by Corsico et al. (2007a) for the hottest known 
GW Vir star RX J2117.1-h3412, and the second being the 
study performed by Corsico et al. (2007c) for the coolest 
member of the class, PG 0122-f200. 

The paper is organized as follows: in the next Section we 
briefly describe our PG1159 evolutionary models. In Sect. [3] 
we derive the stellar mass of PG 1159—035 by means of the 
observed period spacing. In Sect. [J] we derive structural pa- 
rameters of this star by employing the individual observed 
periods. In this section we derive an asteroseismological 
model representative of PG 1159—035 fi i4.3p . and discuss 
aspects such as the mode-trapping properties ( ij4.4p and 
the rates of period changes fi j4.5p of the best-fit model, and 
the asteroseismological distance of PG 1159—035 (i j4.6p . 
Finally, in Sect. [5] we summarize our main results and make 
some concluding remarks. 

2. Evolutionary models and numerical tools 

The pulsation analysis presented in this work relies on a 
new generation of stellar models that take into account the 
complete evolution of PG1159 progenitor stars. Specifically, 
the stellar models were extracted from the evolutionary cal- 
culations recently presented by Althaus et al. (2005), Miller 
Bertolami & Althaus (2006), and Corsico et al. (2006), who 
computed the complete evolution of model star sequences 
with initial masses on the ZAMS in the range 1 — 3.75M0. 
All of the post-AGB evolutionary sequences computed with 
the LPCODE evolutionary code (Althaus et al. 2005) have 
been followed through the very late thermal pulse (VLTP) 
and the resulting born-again episode that give rise to the 
H-deficient, He-, C- and 0-rich composition characteristic 
of PG1159 stars. The masses of the resulting remnants are 
0.530, 0.542, 0.556, 0.565, 0.589, 0.609, 0.664, and 0.741Mo. 
The evolutionary tracks in the log Tcs — log g plane for the 
PG1159 regime are displayed in Fig.[TJ 

The use of these evolutionary tracks constitutes an 
improvement with respect to previous asteroseismologi- 
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Fig. 2. The dipole (left) and quadrupole (right) asymptotic period spacing (AII^) for different stellar masses in terms 
of the effective temperature. The stages before (after) the evolutionary "knee" are depicted with dashed (solid) lines. 
Numbers along each curve denote the stellar masses (in solar units). The plot also shows the location of PG 1159—035 
according to CEA07. 



cal studies. These evolutionary sequences reproduce (1) 
the spread in surface chemical composition observed in 
PG1159 stars, (2) the short born-again times of V4334 Sgr 
(Miller Bertolami ct al. 2006 and Miller Betolami & Althaus 
2007a), (3) the location of the GW Vir instability strip in 
the log Teff — log g plane (Corsico et al. 2006), (4) the ex- 
pansion age of the planetary nebula of RX J2117. 1+3412 
(Corsico et al. 2007ab), and (5) the period spectrum of the 
coolest GW Vir star, PG 0122+200 (Corsico et al. 2007c). 

We computed i — 1 and i = 2 g-mode adiabatic pul- 
sation periods with the same numerical code and methods 
we employed in our previous works (see Corsico & Althaus 
2006 for details). We analyzed about 3000 PG1159 mod- 
els covering a wide range of effective temperatures (5.4 
log(reff) ^ 4.8) and luminosities (0 ^ log(L,/Lo) ^ 4.2), 
and a range of stellar masses (0.530 < M^,/Mq < 0.741). 

3. Mass determination from the observed period 
spacing 

In this section we constrain the stellar mass of 
PG 1159—035 by comparing the asymptotic period spac- 
ing and the average of the computed period spacings with 
the observed period spacing ( ij3.1l and 13.21 respectively). 
These methods take full advantage of the fact that the pe- 
riod spacing of PG1159 pulsators depends primarily on the 
stellar mass, and weakly on the luminosity and the He-rich 
envelope mass fraction (KB94; Corsico & Althaus 2006). 

Most of the published asteroseismological studies on 
PG1159 stars rely on the asymptotic period spacing to infer 
the total mass of GW Vir pulsators, the notable exception 
being the works by KB94 for PG 1159-035, Corsico et al. 
(2007a) for RX J2117.1+3412 and Corsico et al. (2007c) for 
PG 0122+200. 

The precise evolutionary status of PG 1159—035 is not 
known a priori. Thus, we assume two possible cases in our 
following analysis: one situation in which the star is still on 
the rapid contraction phase before reaching its maximum 
effective temperature (i.e. before the evolutionary "knee") 



and the other situation in which the star has just settled 
onto the white dwarf cooling track (i.e. after the evolution- 
ary "knee"). 

3.1. First method: comparison with the asymptotic period 
spacing, An| 

Fig.[2]displays the asymptotic period spacing for £ = 1 (left) 
and £ — 2 (right) modes as a function of the effective tem- 
perature for different stellar masses. Also shown in these di- 
agrams is the location of PG 1159—035, with Toff — 140 + 5 
kK (Werner et al. 1991), and AHg^^ = 21.43 + 0.03 s and 
An^^2 = 12.38 + 0.01 s (CEA07). The asymptotic pe- 
riod spacing is computed as An| = Ho/ \/£{£ + 1), where 
Ho ^ 2Tr'^[J^^{N/r)dr]-\ and N the Bnmt-Vaisala fre- 
quency (e.g. Tassoul et al. 1990). 

From the comparison between the observed AII^ and 
An^ we found a stellar mass between « O.585M0 (if 
the star is located before the evolutionary knee) and « 
O.577M0 (if the star is located after the evolutionary knee), 
irrespective of the £ value; see Table [TJ The inferred range 
of mass is in agreement with the value « O.59M0 de- 
rived by WEA91 and KB94 — and also in agreement with 
the value derived in CEA07 from the KB94 parameteriza- 
tion — on the basis of an asymptotic analysis of the period 
spacing. 

As in our previous works, we must emphasize that the 
derivation of the stellar mass using the asymptotic period 
spacing may not be entirely reliable in pulsating PG1159 
stars (see Althaus et al. 2007). This is because the asymp- 
totic predictions are strictly valid in the limit of very high 
radial order (long periods) and for chemically homogeneous 
stellar models, while PG1159 stars are supposed to be 
chemically stratified and characterized by strong chemical 
gradients built up during the progenitor star life. A more 
realistic approach to infer the stellar mass of PG1159 stars 
is to compare the average of the computed period spacings 
with the observed period spacing. 
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Fig. 3. Same as Fig. [2l but for the average of the computed perfod spacings (AII^). 



Table 1. Stellar mass (in Mq) 
from the period-spacing data. 



of PG 1159-035 derived 





From 


An? 


From 






e = i 


1 = 2 


e = i 


1 = 2 


Before knee 


0.585 


0.585 


0.586 


0.587 


After knee 


0.577 


0.577 


0.561 


0.569 



3.2. Second method: compariso n wit h the average of the 
computed period spacings, AH^ 



The quantity Allf is assessed by averaging the computed 
forward period spacings (Allfc = Ilfc+i — life) in the range 
of the observed periods in PG 1159-035: 390 - 990 s for 
^ = 1 and 350 - 860 s for ^ = 2 (see Tables 4 and 12 of 
CEA07). 

In Fig. [3] we show the run of average of the computed pe- 
riod spacings for £ = 1 (left) and for £ = 2 (right) in terms 
of the effective temperature for all of o ur PG 1159 evolu- 
tionary sequences. Note that the run of Allf depends on 
the range of periods on which the average of the computed 
period spacing is done. Again, the stages before (after) the 
evolutionary knee are depicted with dashed (solid) lines. 
By adopting the effective temperature of PG 1159—035 as 
given by spectroscopy (Teg = 140 ±5 kK) we found a stellar 
mass in the range « 0.586 — O.587M0 (if the star is located 
before the evolutionary knee) and « 0.56 — 0.57Mq (if the 
star is located after the evolutionary knee) ; see Table [1] 
Note that these values are somewhat different to the mass 
derived in Corsico et al. (2006) (« O.558M0) because in 
that paper the authors used a different range of periods 
to compute the average of the period spacing, and a older 
value for observed period spacing value. 

From Tabled] we see that, if the star is evolving before 
the knee, the value of derived from the AII^ is only 
slightly larger (0.25%) than the value inferred from AII^. 
If star is evol ving after the knee, instead, the stellar mass 
derived from AII^ is ^ 3% lower than the mass derived 
from An^. 

In the plot of Fig. [l] the squares correspond to the ap- 
proximate location of the star as predicted by An|, whereas 
the diamonds depict the situation as predicted by AII^ . For 



the case in which the star is evolving before the knee, the lo- 
cation derived from An| coincides with that predicted from 
An^. In this case the star would be located well inside the 
theoretical GW Vir instability strip, but the surface gravity 
would be excessively low (log 5 « 6) as compared with the 
spectroscopically inferred value (log 5 = 7 ± 0.5). Thus, we 
will discard this solution. On the other hand, if the star is 
evolving after the knee, the surface gravity would be some- 
what larger (logg « 7.5) but well within the spectroscopic 
uncertainty, although in this case the nonadiabatic model 
is pulsationally stable against £ = I modes, contradicting 
the observational evidence. 



4. Constraints from the individual observed periods 

In this approach we seek a pulsation model that best 
matches the individual pulsation periods of PG 1159—035. 
The goodness of the match between the theoretical pulsa- 
tion periods (life) and the observed individual periods (IiP) 
is measured by means of a merit function defined as 



X 



n 

■(M^Teff) - - Vmin[(nP -nfe)2]. 



(1) 



1=1 



where n is the total number of the observed peri- 
ods considered. Note that taking the square root of 
we obtain the standard deviation between the ob- 
served and the theoretical periods. The PG 1159 model 
that shows the lowest value of will be adopted as 
the "best-fit model". This approach — which is usually 
called the forward method in asteroseismology — has also 
been used in the context of pulsating PG1159 stars by 
Corsico & Althaus (2006) and Corsico et al. (2007ac). 
We evaluate the function x^(M*,reff) for stellar masses of 
0.530,0.542,0.556,0.565,0.589,0.609,0.664, and O.741M0. 
For the effective temperature we employed a finer grid 
(AToff ~ 100 - 1000 K) which is given by the time step 
adopted in our evolutionary calculations. 

4.1. The WEA91 observed periods 

We first performed period-to-period fits by considering the 
old set ofPG 1159-035 observedm = periods of WEA91. 
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Fig. 4. The inverse of the quahty function corresponding to period fits considering £ = I modes only (left panel), £ — 2 
modes only (centre panel), and £ — 1,2 modes simultaneously (right panel). Lines with different colours correspond to 
sequences with different stellar masses. Dotted (solid) lines correspond to stages before (after) the maximum effective 
temperature of each evolutionary sequence is reached. All of the observed m = periods were extracted from Table 3 of 
WEA91. For clarity, the curves have been artificially shifted upward (with a step of 0.1), except for the lower- most curve 
which corresponds to the O.SSOM© sequence. The vertical grey band displays the PG 1159—035 spectroscopic effective 
temperature and its uncertainties. [Color figure only available in the electronic version of the article] 
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Fig. 5. Same as Fig. SI but for the case of the observed m = periods with identification secure only, extracted from 
Table 3 of WEA91. [Color figure only available in the electronic version of the article] 



They are reproduced in the first and third columns of our 
Table [2] in units of seconds. A letter "s" indicates a mode 
with a secure identification of m and £. The last row shows 
the total number of periods of the corresponding column. 
The results of our period fits are displayed in Figs. [4] and 
O where the quantity (x^)~^ in terms of Tcff for different 
stellar masses is shown. In the interest of clarity, the curves 
are shifted upwards with a step of 0.1 starting from the 
curve corresponding to the O.53OM0 sequence. The spectro- 
scopic Toff and its uncertainties are depicted with a vertical 
grey strip. A peak in the inverse of the merit function with 
{x^)~^ ^ 0-3 will be considered as a good match between 
the theoretical and the observed periods. 

We started our analysis by performing a period fit con- 
sidering the supposed m = 0, £ = 1 modes only (first col- 
umn of Table [21). In this way, the value of £ of the computed 
periods is fixed to be £ = 1 in our period fit procedure. The 
function (x^)~^ is characterized by a rather smooth be- 



haviour except for a few pronounced peaks at stages after 
the turning point in Toff is reached (left panel of Fig. [4]). In 
particular, there exists a primary peak that is seen for the 
first time at high effective temperatures (« 160 000 K) for 
the O.556M0 sequence. The peak gradually shifts to lower 
Toff for the sequences with higher masses, lying at ~ 65 000 
K for M, = O.741M0. The maximum amplitude of the peak 
(i.e., the best period match) is reached for a O.565M0 model 
with Toff ~ 128 000 K, somewhat lower than the spec- 
troscopic effective temperature of PG 1159—035 (Tcff = 
140 000 ± 5 000 K). In other words, if the £ = 1 modes 
were the only modes present in PG 1159—035, our best as- 
teroseismological solution derived from a period fit should 
be a model with Toff = 127 680 K and = 0.565Mq. 
The quality of our period fit is measured by the average 
of the absolute period differences, (511; — (X^ILi l'^nj|)/'^j 
where Slli = Ilf — Ilk, and by the root- mean-square resid- 
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Table 2. Observed periods (11^) of identified m = and 
^ = 1,2 modes for PG 1159-035. 



WEA91 




CEA07 




WEA91 




CEA07 




i = 1 




t = 1 




i = 2 




^ = 2 








390.30 


s 


339.24 


s 










412.01 


s 


351.01 


s 


350.75 




430.04 


s 


432.37 




361.76 


s 


363.39 




450.71 




452.06 


s 






376.03 




469.57 


s 


472.08 




388.07 


s 


387.47 




494.85 


s 


494.85 


s 


398.89 




400.06 


s 


517.18 


s 


517.16 


s 


413.30 


s 


413.14 




538.16 


s 


538.14 


s 


425.03 


s 


425.04 




558.44 


s 


558.14 


s 


438.00 


s 


(437.91) 




581.29 




579.12 




(453.62) 




449.43 




603.04 


s 


603.04 








498.73 


s 


622.60 


s 


622.00 








511.98 


s 


643.41 


s 


643.31 


s 






524.03 




666.22 


s 


668.09 








536.37 




687.71 


s 


687.74 


s 






547.00 




707.92 


s 


709.05 








561.99 


s 


729.50 


s 


729.51 


s 


577.17 


s 


573.69 


s 


753.12 


s 


752.94 








585.26 




773.77 


s 


773.74 




635.92 








790.94 


s 


791.80 




649.67 








817.12 


s 


814.58 








660.46 




840.02 


s 


838.62 
861.72 
883.67 
903.19 


s 
s 


694.88 




672.20 
684.48 
696.83 
709.87 








925.31 


s 


(734.87) 




746.38 


s 






945.01 




779.48 




783.19 








966.98 


s 


812.57 




820.90 


s 






988.13 




833.31 

(929.38) 
982.22 




858.84 




20 


18 


29 


14 


19 8 


26 


7 



Note: parenthesis indicate m = periods whicli are actually 
absent from the power spectrum. Their values are estimated by 
averaging the components m = ±1. 



ual, CT^jj. — V^^l^nTpy/n. For this solution we obtain 
'Wi = 1.32 s and a,^^ = 1.80 s. 

Next, we carried out a period fit taking into account the 
supposed TO = 0, £ = 2 modes only (third column of Table 
[2]). The value of i of the theoretical periods is fixed to be 
£ = 2. The results are depicted in the centre panel of Fig. H) 
The behaviour of (x^)~^ is very different as compared with 
the £ = 1 period fit. Remarkably, (x^)~^ exhibits numerous 
peaks of almost similar amplitude, irrespective of the stellar 
mass. This means that the quadrupole period spectrum of 
the star could not be fitted by a unique PG1159 model, due 
to the existence of numerous and almost equivalent possible 
solutions Q. We conclude that if the 1 = 2 modes were the 

^ This effect can be understood on the basis that the pulsation 
periods of a specific model that evolves towards higher (lower) 
effective temperatures generally decrease (increase) with time. If 
at a given T^s the computed periods of the model show a close fit 
to the observed periods, then the function (x^)~^ reaches a local 
maximum. Later when the model has evolved enough (heated or 
cooled), it is possible that the accumulated period drift nearly 
matches the period spacing between adjacent modes (|Afc| = 1). 
In these circumstances, the periods of the model are able to fit 
the periods of the star again, as a result of which (x^)~^ shows 
other local maxima. 



only present in PG 1159—035, a satisfactory asteroseismo- 
logical solution could not be found. 

Finally we made a period fit using the to = 0, £ = 1 
and to = 0, i = 2 modes simultaneously (first and third 
columns of Table H]) . In this case the value of £ for the the- 
oretical periods is not fixed but instead is obtained as an 
output of our period fit procedure, although the allowed 
values are 1 and 2. The results are displayed in the right 
panel of Fig. |4l The (x^)~^ function shows the same be- 
haviour as in the case of the i = 2 period fit. Thus, we 
are unable to find an asteroseismological model satisfying 
simultaneously both the £ = 1 and £ = 2 sets of observed 
periods. 

Our next step was to repeat the above period fits, but 
this time by considering only the subset of observed periods 
for which the £ and to identification is considered as secure 
and free of ambiguities by WEA91. The periods employed 
are labelled with a letter "s" in Table [2l As a result of this 
selection, the number of observed £ = 1 periods is slightly 
reduced from 20 to 18, but for £ = 2 modes the number is 
substantially reduced from 19 to 8 periods. 

The results of the period fits are shown in Fig. [H To 
begin with, for £ = 1 modes (left panel) the behaviour of 
{x'^)~^ remains unaltered as compared with the situation 
analyzed in the left panel of Fig. 21 due to both period fits 
being essentially the same. For the best-fit model, which is 
again the 0.565Mq model at Toff = 127 680 K, we obtain 
SHi = 1.24 s and a^^, = 1.76 s. For the £ = 2 period fit, on 
the other hand, the (x^)~^ function experiences appreciable 
changes as compared with the situation analyzed before. In 
particular, the numerous peaks present in (x^)~^ exhibit 
larger amplitudes, revealing period matches that are some- 
what better as a result of the smaller number of periods to 
be fitted. 

Finally, the simultaneous fit to ^ = 1 and £ = 2 observed 
periods results again in a (x^)^^ function with a complex 
peaked structure. Nevertheless, it is possible to distinguish 
in this case a prominent peak corresponding to the same 
best-fit model previously found, with Tcff = 127 680 K and 
A/» = O.565A/0, and the period fit in this case is character- 
ized by SHi = 1.25 s and a^^j, = 1.6 s. A detailed compari- 
son of the observed periods with the theoretical periods for 
this period fit is provided in Table [3) 

Note that the assignation of the £ value for the theoret- 
ical periods differs from the identification of the observed 
periods in some cases. 

4.2. The CEA07 observed periods 

In this section we repeat the above analysis but this time 
taking full advantage of the updated and augmented set 
of observed periods of PG 1159-035 reported by CEA07. 
Indeed, these authors have identified 76 new pulsation 
modes, increasing to 198 the total number of pulsation pe- 
riods of PG 1159—035. We reproduce in the second and 
fourth columns of Table [5] the subset of £ = 1 , to = and 
f = 2, TO = periods we use here in our period fits. They 
correspond to Tables 4 and 12 of CEA07. In our Tabled a 
letter "s" means a Confidence Level of 1 or 2. The observa- 
tional uncertainties in the values of the periods are between 
^ 0.002 and ~ 1 s, with an average of ~ 0.3 s. 

^ In the derivation of the quoted errors in the observed peri- 
ods we are assuming that (1) each period corresponds to a real 
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Fig. 6. Same as Fig. 21 Here, all of the observed m = periods (having a Confidence Level of 1, 2, 3, 4, 5, or 6) have 
been extracted from Tables 4 and 12 of CEA07. [Color figure only available in the electronic version of the article] 
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Fig. 7. Same as Fig. [HI but for the case of the observed m — periods with Confidence Level of 1 or 2 only, extracted 
from Tables 4 and 12 of CEA07. [Color figure only available in the electronic version of the article] 



We note that the pulsation periods of PG 1159—035 
are changing with rates up to ^ ±30 ms/yr (Costa & 
Kepler 2007 [CK07]). So, between 1983 and 2002, the pe- 
riods could have experienced variations up to 0.6 s. Thus, 
for modes present in more than one periodogram, the peri- 
ods in Tables 4 and 12 of CEA07 (and also in our Table [l]) 
are the average values. At variance with this, the periods 
included in Table 6 of CEA07 are the values at the epoch 
of the largest amplitude of each mode. Thus, these periods 
correspond to different annual data sets. In our period fit 
procedure below — and also in the analysis of the period 
spacing of Sect. 14.41 — we shall employ the set of the aver- 
age values of periods (Tables 4 and 12 of CEA07) because 
the use of periods from different epochs (Table 6 of CEA07) 
seems inappropriate. 

Regarding the reliability of the detected modes, CEA07 
classified the possible pulsation periods by their confidence 
level (CL), in six different levels of decreasing probability 
from 1 to 6. The six probability levels are listed in Table 2 
of CEA07. 



The results for the i — 1 period fit are displayed in the 
left panel of FiglS] As previously, a quite prominent peak 
in (x^)~^ is found, corresponding to the 0.565Mo model at 
Teff = 127680 K. Remarkably, the period match turns out 
be of an excellent quality, with dUi = 1.03 s and a^-^ — 1.7 
s, in spite of the large number (29) of £ ~ 1 fitted periods. 
The recurrent presence of this solution throughout all of our 
period fits is striking and indicates that this PG1159 model 
could constitute a meaningful asteroseismological solution 
for PG 1159-035. 

The results of the period fit for £ — 2 (centre panel of 
Fig. [S]) are not very different from those of the £ = 2 period 
fits previously discussed. However, in the present case there 
exists a very pronounced peak in (x^)~^ corresponding to 
a O.565M0 model at « 167 000 K, associated with a period 
match characterized by 6Ili — 0.60 s and tr^j^ = 0.33 s. 
Unfortunately, since this period fit relies on the observed 
£ = 2 periods only, which generally have more uncertain 
identifications of £ and m than dipole periodf0, and because 
this solution does not also appear for the dipole modes 



eigenmode of the star, and (2) the indexes k, £, m axe correctly 
assigned. 



^ In particular, the m identification for ^ = 2 is insecure, be- 
cause generally not all the components of the quintuplets are 
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Table 3. Comparison between the observed {£, m) — {1, 0) 
and {£,m) = (2,0) periods [Uf , in units of seconds) 
for PG 1159-035 taken from WEA91 and the theoretical 
{i,m) = (1,0) and {i,m) — (2,0) periods (11^, in units of 
seconds) of the best-fit model with — O.565M0 and 
Teff = 127 680 K. SHi = nf - represents the period 
differences, k the radial orders, and 11^ the rates of period 
change (in units of 10~^^ s/s). 
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— for which we are confident of the i and m values — we 
must discard this model as a meaningful asteroseismological 
solution for PG 1159—035. Other important peak seen in 
(X^)"^ is located at T^s « 100 000 K for AU = O.6O9M0. 
In decreasing order of importance, we found other peaks 
at Teff « 138 000 K (M, = O.565M0), T^ff w 112 500 K 
(M* = O.589M0), Tcff « 83 500 K (M* = 0.664Mq), and 
Teff « 271 000 K (M, = O.74IM0). 

All of the significant peaks described above are ironed 
out when we consider a period fit for the £ = 1 and £ — 2 
periods simultaneously, as can be seen in the right panel 
of the Fig. [51 Again, no clear solution emerges satisfying 
simultaneously both the £ — 1 and £ = 2 sets of observed 
periods, although there is a wide bump just below ~ 140 000 
K for the sequence with Af* = O.565M0. 

Next, we repeated the period fits described above, but 
this time by considering only the most probable modes of 
the data set, i.e. the modes with CL= 1 or 2. Modes with 
CL= 1 have large amplitudes and appear in one or more of 
the periodograms of CEA07, whereas modes with GL— 2 
have lower amplitudes but even above the detection limit, 
and appear in two or more periodograms. The periods of the 
modes with CL= 1 or 2 according CEA07 are labelled with 
a letter "s" in Tabled] Note that in this case the periods 



detected. This leads to large uncertainties in the precise value 
of m for a given mode. 



Table 4. Observed {£,m) = (1,0) periods (nf, in 
units of seconds) and amplitudes (in units of mma) for 
PG 1159-035 taken from CEA07, theoretical {£,m) = 
(1,0) periods (11^ in units of seconds) of the best-fit model 
with = 0.565Mq and Tes = 127680 K. SU, = Hp - Ilk 
represents the period differences, k are the radial orders, ITfc 
are the rates of period change (in units of 10~^^ s/s), and 
Ilfc are the second order rates of period change (in units of 

10-23 S/S2). 
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to be fitted are drastically reduced in number, from 29 to 
14 for ^ = 1 and from 26 to 7 for ^ = 2. 

The results of our analysis are presented in Fig. [71 
Notably, for the £ = 1 period fit we again recover the per- 
sisting solution with Af, = 0.565AfQ and Tcs = 127 680 
K. Remarkably, the quality of this fit for PG 1159—035 is 
much better than of those previously discussed in this pa- 
per for £ ~ 1 modes, with Slli = 0.64 s and cr^n — 0.75 s. 
A comprehensive comparison of the observed periods with 
the theoretical periods for this fit is provided in Table [H 

For the £ = 2 period fit we obtain similar results than 
above (see centre panel of Fig. [7]). In particular, the acute 
peak corresponding to the O.565Af0 model at « 167 000 
K persists, and in the present case it is characterized by 
SHi = 0.71 s and 0-^^ — 0.75 s. The new feature is the 
presence of a pronounced peak in (x^)~^ corresponding to 
a O.664Af0 model at T^g ~ 147000 K, before the turning 
point in Toff is reached. The peak at Teff w 100 000 K for 
Af, — O.6O9Af0, on the other hand, has suffered from a 
notable reduction of its amplitude. 

Again we are unable to found a unambiguous solution 
for the case in which the £ = 1 and £ = 2 sets of observed 
periods are compared simultaneously with £=1,2 theo- 
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Fig. 8. The inverse of the quaUty function corresponding to period fits considering i — 1,2 modes simuhaneously. 
Observed £ — 1 {i = 2) periods are compared with theoretical i = 1 {£ — 2) periods. The left panel corresponds to a 
period fit that includes all of the observed m = periods (having a Confidence Level of 1, 2, 3, 4, 5, or 6) extracted from 
Tables 4 and 12 of CEA07. The right panel corresponds to a period fit in which only the high confidence periods (those 
having CL= 1 or 2) are employed. [Color figure only available in the electronic version of the article] 
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corresponding 



as Fig. [6l but for the case in which the observed m = periods are weighted with the square of the 
amplitude in the period fit procedure. [Color figure only available in the electronic version of the article] 



retical periods. This is evident from the right panel of Fig. 

El 



the other only 21 — lead to similar results lends confidence 
to this asteroseismological solution. 



4.2.1. Fixing the £ value of the observed periods 

Here we perform additional period fits for ^ = 1 and £ = 2 
simultaneously but forcing the £ — 1 observed modes to be 
fitted hy £ — 1 theoretical modes, and the £ = 2 observed 
modes to be fitted hy £ = 2 theoretical modes. In other 
words, we assume that the £ values adopted by CEA07 for 
the observed periods are correct. Fig. [5] displays in the left 
panel the results for a period fit considering all the periods 
(CL= 1, 2, 3, 4, 5, or 6) and in the right panel the results 
for a period fit in which only the CL= 1 or 2 are taken into 
account. Note that in both cases we found again the best- 
fit model with = 0.565Mq and Tes = 127 680 K. The 
fact that both period fits — one involving 55 periods and 



4.2.2. Weighting the period fits 

We have finally performed period fits by using the observed 
periods of CEA07 with CL= 1, 2, 3, 4, 5 and 6, but weight- 
ing each period in the quality function with the corre- 
sponding amplitude. Specifically, the quality function now 
reads: 



x'(M„T;ff) 



E:Umin[(np-nfe)^]A| 



(2) 



where we have adopted the amplitudes Ai of the observed 
modes — extracted from Tables 4 and 12 of CEA07 — as 
the weights. In this way, the period fits are more influenced 
by modes with large amplitudes than by modes with low 
amplitude. In other words, modes of large amplitude must 
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match better the theoretical model than modes of small 
amplitude. Note that if we set Ai = 1 for i = l,---,n, 
i.e., assume the same weight for all of the modes, Eq. ^ 
reduces to Eq. ([1]), which we used in the period fits of the 
previous sections. 

As before, in the period fits the harmonic degree was 
allowed to adopt the values £ = 1 ot £ = 2. The results 
are shown in Fig. [S] Unfortunately, no apparent solution is 
found, neither in the case of ^ = 1 (left panel) or ^ = 2 (cen- 
tre panel) period fits, nor in the case of period fits to ^ = 1 
and £ = 2 modes simultaneously (right panel) . This is some- 
how expected since by using Eq. ^ the period fit is domi- 
nated by only a few large-amplitude modes (those with pe- 
riods 390.30, 452.44, 517.16, 558.14, 791.80, and 883.67 s), 
whereas the remaining modes, which have notably smaller 
amplitudes, barely contribute to the weighting process. As 
a result, multiple and almost equivalent minima in the 
function are obtained, rendering it virtually impossible to 
find an unambiguous asteroseismological solution. 

4.3. The asteroseismological best-fit model 

The results described in the last sections strongly suggest 
the existence of a significant asteroseismological solution 
for PG 1159—035 corresponding to a model with a stellar 
mass of = O.565M0 and T^s = 127 680 K. We arrive 
at this conclusion by considering mostly £ = 1 period fits, 
although this model also provides a good period fit for £ — 1 
and £ — 2 modes simultaneously if we force observed £ — 1 
periods to be compared with theoretical £ = 1 periods, and 
observed £ = 2 periods to be compared with theoretical 
£ = 2 periods (Sect. Bmj) . Also, we found a good period fit 
for this model when we consider the observed periods with 
secure identification of WEA91, allowing the periods to be 
^ = 1 or = 2 without constraining any to be one or the 
other (Table [3]). This model will be adopted as the "best-fit 
model" representative of PG 1159-03O 

The quality of our ^ = 1 period fits for PG 1159-035 
using the periods of CEA07 is characterized by Slli = 1.02 
s (when we consider modes with CL= 1, 2, 3, 4, 5, or 6; 
see Fig.|6]) and SHi — 0.62 s (if we take into account modes 
with CL= 1 or 2 only; see Fig. [7]). These period fits are of 
better quality tha n tho se obtained by K B94 and Corsico 
& Althaus (2006) {Wi = 1.19 s and JIlI = 1.79 s, respec- 
tively) for PG 1159-035. The £ = 2 periods of our best-fit 
model do not match the observed quadrupole period spec- 
trum as accurately as the dipole periods do. In spite of this 
fact, the best-fit model yields an average of the computed 
period spacings for f = 2 of An^=2 = 12.61 s, in good agree- 
ment with the observed value of An2^2 — 12.38 ± 0.01 s. 
This property renders consistency to the derived best-fit 
model. 

As in our previous asteroseismological analysis of GW 
Vir stars, we are able to get a PG1159 model that repro- 
duces the period spectrum observed in the star under study 
without tuning ad hoc the value of structural parameters 
such as the thickness of the outer envelope, the surface 



It is important to stress that period fits using only £ = 2 
modes show interesting possible solutions, in particular the 
O.565M0 model at ^ 167000 K described in the last section. 
However, we must discard these solutions because they exclude 
the £ — 1 modes (which have more secure identifications than 
the £ — 2 modes). 



chemical abundances or the shape of the core chemical pro- 
file which, instead, are kept fixed at the values predicted by 
the evolutionary computations. 

Table |4] shows a detailed comparison between the 29 
£ — 1 CEA07 observed periods (Hp) and the theoretical 
£ — 1 periods (11^) of our best-fit model. Note that if we 
consider all the 29 periods (CL= 1, 2, 3, 4, 5 or 6), we have 
l^n^l < 1 s for 17 modes (59%), 1 s < l^H^I < 2 s for 7 
modes (24%), 2 s < l^n^l < 3 s for 4 modes (14%), and 
l^IIil > 3 s for 1 mode (3.5%). However, if we adopt the 
14 periods with the highest level of probability only (CL= 
1 or 2), we obtain \6Ili\ < 1 s for 11 modes (« 80%) and 
1 s < I^HjI < 2 s for 3 modes (« 20%). 

The fifth column of Table |4] presents the radial order 
k for the modes of our best-fit model. The k values are 
in complete agreement with those of the best-fit model of 
KB94, and also in agreement (to within Afc = ±2) with 
the values proposed by CEA07. As a reference, we show 
the radial order of the mode with the largest amplitude 
(n ~ 517 s), according to the different authors, in the last 
row of Table El 

In Table [5] we show the main features of our best-fit 
model and also the parameters of PG 1159—035 extracted 
from other published studies. Specifically, the second col- 
umn corresponds to spectroscopic results, the third col- 
umn presents results from the asteroseismological analysis 
of KB94, the fourth column corresponds to the pulsation 
resuhs of CEA07 and Costa & Kepler (2007) (CK07), and 
the fifth column shows the characteristics of the asteroseis- 
mological model of this work. The quoted uncertainties in 
the stellar mass and the effective temperature of our best 
fit model (ctm, and aT^„) resulting from the period fit pro- 
cedure are calculated according to the following expression, 
derived by Zhang et al. (1986): 



where So — x^{M^J,T^^) is the absolute minimum of 
which is reached at (M^,TiL) corresponding to the best-fit 
model, and S the value of x when we change the parameter 
i (in our case, M* or Tes) by an amount di keeping fixed the 
other parameter. The quantity di can be evaluated as the 
minimum step in the grid of the parameter i. Specifically, 
we use the £ = 1 period fit shown in the left panel of Fig. 
[3 that is, by considering the 14 dipole periods with CL= 
1 or 2 measured by CEA07. We have dT^f = ATes ~ 1000 
K and dM, = AM, in the range 0.009 -°O.O24M0. 

Our best-fit model has a stellar mass of = 0.565Mo, 
which is ^ 1.9% smaller than the value derived from the 
average of the computed period spacing, M, ~ 0.576Mo, 
and ~ 3% lower than that inferred from the asymptotic 
period spacing, ~ O.581M0 (see fj3]). On the other 
hand, CEA07 have inferred a value of the stellar mass of 
PG 1159—035 by using an interpolation formula to the pe- 
riod spacing derived by Kawaler & Bradley (1994) on the 
basis of a large grid of artificial PG1159 models in the lumi- 
nosity range 1.6 ^ log(i,/L0) <: 3.0. These authors obtain 
a rather high value of 0.59 ± 0.02Mq, in line with the early 
determinations of WEA91 and KB94 (0.586 ±0.003Mo and 
O.59±O.O1M0, respectively), but in disagreement — about 
4.5% larger — with the mass of our best-fit model. 

The value of our best-fit model is somewhat higher 
than the spectroscopic mass of O.54M0 derived by Miller 
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Quantity 


Spectroscopy 


Asteroseismology 
(KB94) 


Pulsations 
(CEA07, CK07) 


Asteroseismology 
(this work) 


Teft [kK] 


140 ± 5^^^ 


136 


— 


128l|^ 


M, [Mq] 


0.54 iO.l'*'' 


0.59 ±0.01 


0.59 ±0.02 


0. 565^0 


loe fcm/s^l 


7.0 ± 0.5'^' 


7.38 ± 01 




7.42+»?i 

— 0.12 


los(L./Ln) 


^- ' ^-0.55 


2.29 ± 0.05 




2.15 ± 0.08 


loK(ii./-Rn) 




-1.60 ± 0.09 




-1.62+°"g 


Menv [Mq] 





0.002 


— 


0.017 






0.60 - 0.65 


0.83 ±0.05 


0.55 


C/He, 0/He(*' 


1.5, 0.5^"^ 






0.81,0.56 


My [magj 


— 


— 


— 


6.975 ± 0.4 


Mboi [mag] 








—0.625 ± 0.2 


Ay [mag] 








Q Qf;+0.003 
^■^"-0.004 


d [pc] 


onn+«oo(a) 

O00_4oo 


440 ± 40 




363+61 


TT [mas] 


1 nr + 1.25 


o o'7+0.23 
^•^ ' -0.18 




9 7r+0.56 


R,/R, [10-"s-^] 
f/T [10~"s~^] 






-0.89 ±0.2 
-1.84 ±0.04 


-0.009 
0.0004 


k for n = 517 s 




22 


20 


22 



Note: (*) Abundances by mass 
References: (a) Werner et al. (1991); (b) Miller Bertolami & Althaus (2006); (c) Werner & Herwig (2006). 
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Fig. 10. Distribution of the forward period spacing All 
for the case of the observed periods according to Table 6 
— panel (a) — and according to Table 4 — panel (b) — of 
CEA07. 

Bertolami & Althaus (2006) (see also Dreizler & Heber 
1998) for PG 1159—035. Note that a discrepancy between 
the asteroseismological and the spectroscopic values of 
M^, is generally encountered among PG1159 pulsators (see 
Corsico et al. 2006, 2007ac and Althaus et al. 2007). Until 
now, the asteroseismological mass of PG 1159—035 has 
been about 9% larger (AAf* w 0.05Mq) than the spec- 
troscopic mass. In light of the best-fit model derived in this 
paper, this discrepancy is reduced to less than about 5% 
(AM* « 0.025Mq). 

The location of the best-fit model in the log TeS — log g 
diagram is displayed in Fig. [TJ Note that the model is char- 
acterized by a Tcfi somewhat lower than the spectroscopic 
value. Also, the model lies certainly slightly outside the 
£ = 1 GW Vir instability domain, at odds with the obser- 
vational evidence. 



4.4. Period spacing and mode trapping 

The mode-trapping properties of PG1159 stars have been 
discussed at length by KB94 and Corsico & Althaus (2006); 
we refer the reader to those works for details. We will con- 
sider the observed i = 1 forward period spacing (AH = 
Ilfc+i — Ilfc) taking advantage of the uninterrupted string 
of 29 dipole periods with consecutive radial orders from 
fc = 14 to fc = 42, according CEA07. 

In panel (a) of Fig. [TO] we show the AH vs 11 diagram 
for the observed periods according to Table 6 of CEA07. 
The uncertainties in the observed period spacings are be- 
tween 0.01 and 1.2 s, with an average value of 0.4 s. As 
we mentioned at the beginning of Sect. 14. 2[ the periods of 
PG 1159—035 are changing with time. So, the values in- 
cluded in Table 6 of CEA07 are the periods corresponding 
to the epoch of the largest amplitude of each mode. This 
diagram is used by those authors to identify five modes 
trapped in the He-rich envelope of PG 1159—035. They ar- 
rive at this conclusion under the assumption that trapped 
modes appear as points of minimum in this diagram. By as- 
suming the presence of a single chemical interface, CEA07 
derive Tc/R* ~ 0.83, where Tc is the radial coordinate 
of the O/C/He chemical interface. In our best-fit model 
this chemical interface is located at markedly deeper lay- 
ers (rc/i?* ~ 0.55), in agreement with the models of KB94 
(rc/i?* « 0.60-0.65). 

Note that some periods used to construct the diagram 
of panel (a) in Fig. [TO] (extracted from Table 6 of CEA07) 
slightly differ from the periods included in Table 4 of 
CEA07 (and in our Tables [5] and H]) because, as stated be- 
fore, they are the average values of periods on several an- 
nual data sets. The slight differences in some of the periods 
are translated into appreciable differences in some values 
of the period spacing, as can be seen in panel (b) of Fig. 
ITOl where we plot the same diagram but on the basis of the 
periods given in Table 4 of CEA07. Since the period fits 
performed in the present work are based on these periods, 
we shall concentrate on the diagram of panel (b) of Fig. [TO] 
It can be compared with panel (a) of Fig. [TT] in which the 
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Fig. 11. Panel (a): distribution of AH for the best-fit 
model, the primary minima being marked with circles. 
Panel (b): the computed distribution of the kinetic energy. 
Panel (c): the theoretical values of the relative rates of pe- 
riod change. Modes trapped in the outer envelope model 
are labelled with their radial order k. 



An vs n diagram for the best-fit model is displayecH. Note 
that, in spite of the excellent match between the observed 
periods and the theoretical periods of the best-fit model 
(Table m, the pertinent period spacing distribution looks 
somewhat different. The disagreement has its origin in the 
fact that small differences in period — because the period 
match is not perfect — are translated into large differences 
in period spacing. Besides the disagreement between indi- 
vidual values of An, we note that generally the amplitude 
of the observed period spacing are remarkably larger than 
that of the best-fit model. This suggests that the chemical 
interfaces present in PG 1159—035 could be considerably 
steeper than predicted by our PG1159 modeling. 

Notwithstanding the mentioned differences between the 
period spacing distributions in PG 1159—035 and in our 
best-fit model, they share a common feature: both ex- 
hibit primary and secondary minima, characteristic of mode 
trapping due to the presence of more than one chemical in- 
terfaces. In fact, our best-fit model has two chemical tran- 
sition regions: the inner interface of 0/C and the more ex- 
ternal interface of 0/C/He. As we have already demon- 
strated (see, e.g., Corsico & Althaus 2006 and Corsico et 



^ The uncertainties in Tcs and M* of the best-fit model (see 
Table [S| lead to uncertainties in the computed AH, with an 
average value of ~ 1 s. In the quoted errors of T^a and M, 
of the best-fit model we are neglecting possible uncertainties 
in the stellar evolution and pulsation calculations. Thus, the 
quoted error in AH must be taken as lower limit. 



el. 2007ac), for periods below « 650 — 700 s the mode- 
trapping features of our PG1159 models are induced mostly 
by the chemical gradient at the 0/C/ne interface, the 0/C 
interface being more relevant for longer periods. We con- 
clude that the presence of primary and secondary minima 
in the An^ pattern of PG 1159—035 is an indication that 
the star could be characterized by more than one chemical 
interfaces, possibly twcl^. 

Regarding the identification of trapped modes by means 
of the An vs n diagram — which constitutes the main ob- 
servational tool of mode trapping in PG1159 stars — we 
emphasize that generally a trapped mode with radial order 
k is associated with a minimum An which corresponds to a 
radial order fc± 1, and in some cases to k. This is illustrated 
in panel (b) of Fig. II 11 in which we show the kinetic energy 
distribution i?kiiJj, of our best-fit model. Local minima in 
the -Bkin distribution correspond to trapped modes, because 
these modes propagate mainly in the low density regions 
of the outer ne-rich envelope. Note that, for instance, the 
minimum of An with radial order fc = 23 is related to a 
minimum in i?kin of the trapped mode with k — 24. Modes 
trapped in the envelope are also characterized by minima in 
n/n, because they should be more strongly sensitive to the 
effects of the surface contraction — that induces a secular 
decrease of the periods — than untrapped modes — which 
are more affected by cooling. This effect is well illustrated 
in panel (c) of Fig. [TlJ 

4.5. Period changes and the cooling rate of PG 1159— 035 

In Table |4] we show the rate of period change for the 
modes of our best-fit model (sixth column). Our calcula- 
tions predict all of the pulsation periods to increase with 
time (n > 0) , in accordance with the decrease of the Brunt- 
Vaisala frequency in the core of the model induced by cool- 
ing. This property is shared by all of the PG1159 models 
located below the lower thin dotted line in Fig.[TJ Note that 
at the effective temperature of our best-fit model, cooling 
possibly has the largest effect on II, but gravitational con- 
traction — which should result in a decrease of periods with 
time — is not negligible and could affect the n values, in 
particular for the case of modes trapped in the envelope 
(see Fig.HH). 

Our best fit model does not reproduce the recent mea- 
surements of CK07, which indicate that the pulsation 
modes of PG 1159—035 have positive and negative values 
of n. Our PG1159 models having a mix of positive and 
negative n values are located between the two thin dotted 
lines in Fig.[T] Remarkably, the measurements of CK07 are 
consistent with the spectroscopic location of PG 1159—035 
on the log Toff — \ogg diagram. Thus, to be consistent with 
the findings of CK07 our best-fit model for PG 1159-035 
should be located at earlier evolutionary stages. 

CK07 report values of the rate of period change in 
PG 1159—035 that are generally more than one order 
of magnitude larger than the theoretical II values char- 
acterizing our best-fit model. In particular, they report 
ri = {+18.2 ± 0.8) X 10"" s/s for the large amplitude 



^ This conclusion is valid even if we consider the observational 
uncertainties in AH. 

^ The usual normalization of the radial eigenfunction (j/i = 
^r/r — 1 ai the stellar surface) is assumed. 
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m — 0, £ — I mode with period 517.1 s, excessively larger 
than the value predicted by our computations for that mode 
(n = 1.40 X 10~" s/s). In addition, CK07 have been able 
to measure the second order rates of period change (II) for 
eight modes in PG 1159—035. Again, the measured values 
(|II| ~ 0.01 — 9 X 10~^^ ss~^) are strikingly larger than our 
theoretical predictions (|n| ^ 0.1 — 5 x 10"^'^ ss~^). 

Finally, by using the H/H values measured for the 
(m = 0) 517.1 s and the (m = +1) 516.0 s modes, CK07 
roughly estimate for PG 1159—035 a contraction rate of 
= -8.9 ± 0.2 x 10"i2g-i and a cooling rate of 
f/T = -1.84±0.04 X 10""s"\ orders of magnitude larger 
than our theoretical predictions for the best-fit model, of 
= -9.0 X lO'^'^s'^ and f/T = 4.0 x lO'^^s'^ 

4.6. The asteroseismological distance and parallax of 
PG 1159- 035 

We employ the luminosity of our best-fit model to infer the 
seismic distance of PG 1159—035 from the Earth. First, 
we compute the bolometric magnitude from the luminos- 
ity of the best-fit model by means of Mboi = Mq^oI — 
2.51og(L*/Lo), with Moboi = 4.75 (Allen 1973). Next, we 
transform the bolometric magnitude into the absolute mag- 
nitude. My = Mboi — BC, where BC is the bolometric cor- 
rection. For PG 1159-035 we adopted BC = -7.6 ± 0.2 
from Werner et al. (1991). We account for the interstel- 
lar absorption, Ay, using the interstellar extinction model 
developed by Chen et al. (1998). We compute the seismic 
distance d according to the well-known relation: logd = 
i [mv — Mv -I- 5 — Av(<i)] where the apparent magnitude is 
rriv = 14.84. The interstellar absorption Av(d) varies non 
linearly with the distance and also depends on the Galactic 
latitude (6) and longitude {£). For the equatorial coordi- 
nates of PG 1159-035 (Epoch B2000.00, a = 12^ 1™ 46M, 
6 = —3° 45' 39") the corresponding Galactic coordinates 
are b = -56° 51' 55".44 and £ = 279° 49' 10". 2. We solve 
for d and Ay iteratively and obtain a distance d 360 pc 
and an interstellar extinction Ay ^ 0.06. Note that our 
distance is a factor of more than two smaller than the es- 
timation of 8001400 of Werner et al. (1991) and with its 
accuracy substantially improved. Finally, our calculations 
predict a parallax of tt ~ 2.8 mas. 



5. Discussion and conclusions 

In this paper we carried out a comprehensive asteroseis- 
mological study of PG 1159—035, a g-mode pulsator that 
defines the class of pulsating PG1159 stars — the GW Vir 
variables — and the PG 1159 spectral class of hydrogen- 
deficient stars. Our analysis is based on the full PG1159 
evolutionary models presented in Althaus et al. (2005), 
Miller Bertolami & Althaus (2006), Corsico et al. (2006) 
and Corsico et al. (2007c). These models represent a solid 
basis to analyze the evolutionary and pulsational status of 
GW Vir stars like PG 1159-035. 

We first took advantage of the strong dependence of 
the period spacing of variable PG1159 stars on the stellar 
mass, and derived a value of the mass in the range 0.577 — 
0.585Mq by comparing the observed mean period spacing 
An*-' with the asymptotic period spacing of our models. 
We also compared AII*^ with the computed period spacing 



averaged over the period range observed in PG 1159—035, 
and derived a value in the range 0.561 — 587M0. Note that 
in both derivations of the stellar mass we made use of the 
spectroscopic constraint that the effective temperature of 
the star should be ~ 140 kK. 

Next, we adopted a less conservative approach in which 
the individual observed pulsation periods alone — i.e., ig- 
noring "external" constraints such as the spectroscopic val- 
ues of the surface gravity and effective temperature — natu- 
rally lead to an "asteroseismological" PG1159 model that is 
assumed to be representative of the target star. Specifically, 
the method consists in looking for the model that best re- 
produces the observed periods. The period fits were made 
on a grid of PG1159 models with a quite fine resolution in 
effective temperature although admittedly coarse in stel- 
lar mass. We employed the period data from WEA91 and 
also the more modern — and enlarged — period data of 
CEA07. We considered period fits taking into account ob- 
served £ = 1 modes only, observed £ — 2 modes only, and 
observed £ = 1 and £ — 2 modes simultaneously. We found 
a persisting and significant asteroseismological solution for 
PG 1159—035 corresponding to a model with a stellar mass 
of = 0.565Mq and T^s = 127 680 K whose £ = 1 peri- 
ods match the observed periods identified as £ — 1 with an 
average of the period differences (observed versus theoret- 
ical) of only 0.64 s (for 14 £ — 1 periods fitted) and 1.03 s 
(for 29 ^ = 1 periods fitted). Also, for this "best-fit model" 
we obtain a good period match when we consider a period 
fit for dipole and quadrupole modes simultaneously if we 
constrain the £ value of the input periods from the outset. 

The main results of this work are summarized in Fig. [1] 
and TableO where the best-fit model solution as well as the 
inferences from the period spacing analysis are displayed 
together with the spectroscopic inferences. Note that the 
best-fit model is characterized by a stellar mass similar to 
the mass value predicted by the period spacing approach, 
although with a lower effective temperature (slightly out- 
side the 1 a error bar). However, at this location, and 
according to the nonadiabatic pulsation analysis of the 
PG1159 sequences employed here (Corsico et al. 2006) the 
model should be pulsationally stable against £ = 1 modes 
and the rates of period changes all positive, in disagree- 
ment with the observations. In this sense, we must warn 
the reader that the location of the theoretical blue edge of 
the GW Vir instability strip shown in Fig. [1] is not to be 
taken at face value. This is because changes in the surface 
abundances predicted by stellar modeling are expected to 
shift the predicted blue edge location. For instance, Quirion 
et al. (2007) have found that increasing the carbon and oxy- 
gen abundance in the outer layers shifts the blue edge to 
larger effective temperatures. In particular, by varying the 
oxygen abundance from Xq ~ 0.20 to 0.40 (with He being 
replaced by O) and keeping fixed the carbon abundance at 
Xc ~ 0.40, the blue edge of the GW Vir instability strip 
can be pushed to higher effective temperatures by about 
10 000 K. 

We can independently constrain the location of 
PG 1159—035 in the logTeff — logg plane by examining 
the trend of secular period changes. In fact, CK07 have 
reported the period change of 27 pulsation modes. They 
found that some modes show an increase in their periods, 
while other exhibit a clear decrease. This could be refiecting 
that PG 1159—035 is a pre- white dwarf that still has signif- 
icant ongoing contraction. To elaborate this point further. 
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we mark in Fig. [T] with thin dotted hnes the Umits of the 
region for which our PG1159 sequences predict the simulta- 
neous presence of modes with positive and negative fl[3. To 
be consistent with the findings of CK07 our best-fit model 
for PG 1159—035 should be located at earlier evolutionary 
stages, as suggested by spectroscopy. In fact, our predicted 
best-fit model falls in the region of the log Toff — \ogg di- 
agram where all of the pulsational modes are expected to 
exhibit a positive value of the period change. 

Finally, we note that generally the rates of period 
change measured by CK07 are about one order of mag- 
nitude larger than the values predicted by the best-fit 
model. This is particularly true for the large amplitude 
mode at 517.1 s, for which CK07 report a value n = 
(-f 18.2±0.8) X 10^^^ s/s, at odds with our theoretical value 
n = 1.40 X 10~^^ s/s — the same problem occurs with 
all other previous calculations, including those of Kawaler 
(1986). Also, our best-fit model for PG 1159-035 fails to 
reproduce the large values of the second order rates of pe- 
riod changes observed in the star, and the contraction and 
cooling rates estimated by CK07. Those authors do point 
out that these estimates need corroboration. 

The results of the period-fit procedure carried out in 
this work suggest that the asteroseismological mass of 
PG 1159-035 (- 0.565Mq) could be ~ 4% lower than 
thought hitherto O.59M0; see WEA91, KB94 and more 
recently CEA07) and in closer agreement with the spec- 
troscopic mass of O.54M0 as derived by Miller Bertolami 
& Ahhaus (2006) (see also Dreizler & Heber 1998). This 
suggests that a reasonable consistency between the aster- 
oseismological mass and the spectroscopic mass should be 
expected when the same evolutionary tracks are used for 
both derivations of M*, as we have done in Corsico et 
al. (2007c) for PG 0122-f 200 and we do in the present 
work for PG 1159—035. The exception to this assertion 
is RX J2117. 1-1-3412, for which we found an asteroseis- 
mological mass about 25% lower than the spectroscopic 
value (Corsico et al. 2007ac), employing the same stellar 
evolution and pulsation modeling as in the present work. 
As we suggested in that paper, the discrepancy in mass 
could be due to large errors in the spectroscopic determina- 
tion of logg and Tcs for RX J2117. 1-1-3412. Uncertainties 
in the location of the evolutionary tracks in the HR and 
logTefi — logg diagrams due to the modeling of PG1159 
stars and their precursors can be discarded in view of 
the recent work by Miller Bertolami & Althaus (2007b). 
They conclude that the tracks of PG1159 stars of Miller 
Bertolami & Althaus (2006) are robust enough as to be 
used for spectroscopical mass determinations of PG 1159- 
type stars. 

In closing, in this paper we have been able to find a 
PG1159 model that reproduces the observed pulsation pe- 
riods of PG 1159—035 without invoking any artificial ad- 
justments of the structural parameters (the 0/C chemical 
profile, the surface He abundance, or the thickness of the 
He-rich envelope) which, instead, are kept fixed at the val- 
ues predicted by our evolutionary computations. In partic- 
ular, our PG1159 models are characterized by thick He-rich 
envelopes. 



Above the upper bound of this region, all of the modes are 
expected to have a decrease in their periods. 



However, our best-fit model is unable to explain a num- 
ber of important observed properties of PG 1159—035, such 
as: 

— the nature itself of PG 1159—035 as a variable star, 
because the best-fit model lies outside the theoretical 
dipole GW Vir instability domain, i.e., the nonadiabatic 
treatment indicates pulsational stability of the model 
against £ = I (^-modes, 

— the larger range of the period spacings of PG 1159—035, 
as compared with those of the best-fit model, 

— the mixture of positive and negative rates of period 
change measured in PG 1159—035, because the best- 
fit model has all the theoretical rates of period changes 
positive, 

— the large magnitude of the rates of period change de- 
tected in PG 1159—035, because they are an order of 
magnitude larger than the theoretical values of the best- 
fit model. 

In view of the above shortcomings, we must investigate 
if PG 1159—035 could harbor a thinner helium-rich enve- 
lope than predicted by our evolutionary models, a possibil- 
ity sustained by the fact that PG 1159 and born-again stars 
are observed to suffer from appreciable mass loss. Further 
asteroseismological analysis on the basis of PG1159 evo- 
lutionary models characterized by thin He-rich envelopes, 
and additional observational campaigns will be needed to 
make definite conclusions about the internal structure and 
evolutionary status of PG 1159—035. 
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